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The penetration of grid integrated distributed generation (DG) in the present 
decade, has benefited rural communities, the environment, and the power 
sector. These renewable power sources based DGs could eliminate the need 
of extensive transmission networks, especially in remote areas, reduce 
emissions and improve power supply reliability. A significant drawback of 
grid integrated DG systems is the islanding of DG units, which puts workers' 
safety at risk and raises the possibility of damaging electrical infrastructure. 
Therefore, islanding detection techniques are used to reduce the danger 
associated with islanded functioning of DG units. Fast detection, small non 
detection area and less power quality disturbance are the major requirements 
of any islanding detection method. To address this issue of islanding, 
researchers have proposed various islanding detection strategies. This paper 
compares various q-axis controller-based islanding identification approaches: 
sub-harmonic perturbation (SHP), complementary reactive power 
perturbation (CRPP), and even harmonic perturbation (EHP). In all three 
proposed methodologies, the perturbations introduced result in frequency 
deviations surpassing the predefined threshold values. But the time of 
islanding detection is least in the CRPP approach. CRPP can also drift the 
total harmonic distortion (THD) beyond the corresponding threshold in an 
appreciable way. The performance of these (Islanding detection methods) 
IDMs is evaluated through simulations using MATLAB-Simulink on a PV 
fed DG. The efficacy of the comparative analysis is ensured with necessary 
waveforms. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Praveen Raj Rajaswamy Sarojam 


Department of Electrical and Electronics Engineering 
Mar Baselios College of Engineering and Technology 


Nalanchira, Kerala, India 


Email: india.praveenraj @ gmail.com 


1. INTRODUCTION 


Distributed generation (DG) is essential for powering remote areas distant from centralized power 
stations, like hydroelectric or thermal plants, which pose economic challenges. The continuous evolution of 
DGs, driven by advancements in power electronics, has further boosted their popularity [1], [2]. The significant 
reduction in transmission losses associated with renewable energy sources (RES)-fed DGs is an added 
advantage. Integrating these DGs with the grid is crucial and beneficial [3], [4]. However, before achieving 
synchronization, several challenges need to be resolved. One critical concern is the issue of islanding, which 
can lead to serious problems and even fatalities. To address this, the IEEE 1547 standard mandates that DGs 
must disconnect from the electrical power system (EPS) within a short 2-second timeframe if an islanding 
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situation occurs. According to IEEE1547 standard, islanding occurs when a DG keeps supplying power to an 
area despite disconnection from the bulk EPS [5]. Both voluntary and involuntary islanding states can occur in 
a power system. 

In grid synchronized DG (GSDG), when the local load exceeds power generation, the grid shares the 
deficit power. Conversely, when DG produces more power than needed, the surplus is fed back to the grid. 
This makes GSDG operation necessary and beneficial. During loss of mains (islanding), the DG solely supplies 
the load. Researchers have formulated several remote as well as local islanding detection methods (IDMs). The 
local IDMs can further be categorized into passive, active and hybrid IDMs. The communication between 
utility grid and DG is employed in remote IDMs, to monitor and identify unintended islanding. This method is 
advantageous as it is not affected by local power balance, causes no grid disruptions, and can be customized 
for different grid conditions, eliminating non-detection zone (NDZ) issues. Some examples of remote IDM 
include impedance insertion method [6], power line signal scheme, signal produced by disconnect [7], and 
transfer trip scheme [8]. These IDMs theoretically have no NDZ but can be quite expensive. Passive IDM 
employs various electrical parameters at point of common coupling (PCC), such as voltage variation beyond 
thresholds [9], frequency variation beyond thresholds [10], voltage total harmonic distortion (THD) going 
above the desired limits, current THD going above the desired limits [11], rate of change of reactive power 
[12], harmonic signatures [13] in order to recognize an islanding scenario. As the grid's power sharing 
approaches zero, passive islanding detection methods (IDMs) encounter challenges in discerning islanding 
events due to the diminishing significance of deviations in PCC parameters. This limitation encouraged the 
development of active IDMs. Active IDMs purposefully cause the grid and local load at PCC to be perturbed 
by DG inverters. The low impedance of the DG absorbs these disturbances while it is synchronized with the 
grid. However, during islanding, these injected disturbances result in the divergence of one or more PCC 
parameters from their predefined thresholds, resulting in the relay tripping. High frequency signals injection 
techniques [14], d axis current injection [15], reactive power injection, sandia Frequency [16], are a few 
examples of active IDMs. 

Active IDMs can have unfavorable impact on power quality (PQ) even though they yield less NDZ 
[17]. The integration of passive and active islanding detection methods (IDMs) gives rise to the development of 
Hybrid IDMs. D axis current injection combined with rate of change of frequency (ROCOF) [18], active power 
change as the active IDM with rate of change of voltage as passive IDM [19], Sandia frequency shift (SFS), and 
rate of change of frequency (ROCOF) [20] are some of the hybrid IDMs presented in the literature. Hybrid 
methods enable the utilization of the advantages offered by both approaches, while the limitations inherent to 
each method are not avoidable. 

This paper mainly deals with the comparison of the q-axis perturbation-based IDMs. The d-axis 
disturbances are variations in the active power component of the electrical variables. They are primarily 
associated with changes in the real power component and are more relevant for power system control, stability, 
and load management. The q-axis disturbances cause variations in the electrical variables such as voltage 
magnitude and frequency. q-axis disturbances are highly effective in detecting islanding events, as they are 
less likely to occur during normal grid-connected operation. Different q axis perturbation-based islanding 
identification schemes are summarized in Table 1. 


Table 1. Comparison of various performance parameters of Q axis perturbation-based IDMs 


Ref Non-detection zone (NDZ) Fault detection zone (FDZ) _ Identification time Power quality (PQ) 
[21] Small NDZ NO FDZ 0.17s Affects PQ 
[22] Very small, but non-zero NDZ Asymmetrical fault leads to 0.433 s Affects PQ 
false detection 
[23] No NDZ FDZ to be analyzed 0.1785s Affects PQ 
[24] Very small - 0.25s Affects PQ 
[25] No NDZ NO FDZ 0.12s Affects PQ Slightly 
[26] No NDZ Exists FDZ when grid 0.18s Affects PQ but within limits 
frequency exceeds +0.5Hz 
[27] No NDZ - 0.15s Affects PQ 
[28] NDZ exists - 0.16s Larger PQ issues when NDZ is reduced 


An effective islanding detection scheme should satisfy the following criteria: i) Identify islanding 
occurrences across a diverse range of operational scenarios; 11) Discriminate between islanding and non- 
islanding events; iii) Quick detection time (< 2 seconds) is desirable; iv) Compatible with other existing 
DGs; and v) Simple and cheap in implementation. Therefore, a simple and robust islanding detection method 
with minimum detection time is inevitable in GSDGs. Three active methods, based on perturbations through q 
axis current, are presented and compared in this paper. As the reactive power is a deciding factor of frequency, 
the response to detect islanding state could be ROCOF. 
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The major contributions of this work are: i) Comparative analysis of NDZ and FDZ is performed; ii) 
Comparative analysis regarding THDs during grid integrated mode is carried out; iii) Comparative analysis of 
THDs during islanded mode is also presented; and iv) Detection time for different load quality factors (QF) =5 
is presented for the three methods. The paper is organized as: i) The three active islanding detection methods 
along with its response is presented in section 2; 11) Response of the passive methods presented in section 3; 
iii) Results and waveforms necessary for the analysis are depicted in section 4; iv) Comparison of the three 
approaches is depicted in section 5; and v) conclusion in section 6. 


2. THE PROPOSED SCHEMES FOR ISLANDING DETECTION 

In the proposed work three q axis perturbation based active islanding detection methods are presented, 
analyzed and are compared in various perspectives. As q axis disturbances do not pose any stability issues to 
the DG or utility grid, whereas the d axis perturbations can cause stability problems. The proposed system in 
the block schematic is presented in Figure 1. 


PV Module Boost Converter Inverter Filter 


DATA PROCESSING 
p P* 


Figure 1. Block diagram of the proposed system with three q axis control methods 


2.1. Active islanding detection scheme -I (Sub-harmonic perturbation (SHP)) 

Sub harmonic perturbation (SHP) is injected in the q axis current in the d-q frame. So, the perturbation 
fed into the grid at PCC is in such a level that the grid power quality (PQ) is not deteriorated beyond the 
threshold level. In response to the modifications in the q axis current, DG frequency signatures change. When 
the DG is integrated to bulk EPS, voltage at PCC, Vpcc is same as grid voltage, Veriq and is presented in (1). 
The frequency: 


Voce = Verid (1) 


the d axis and q axis current equations of DG inverter is given by (2), 


ae faa (2) 
Iginv 
The current equations are altered as follows when SHP is injected as in (3) and (4). 
ig idinv 
F 3 
f ‘|= f Iginv + Liver 0) 
igpert = = iginvSin(@pertt) (4) 


iginy is DG inverter current. After the Park’s transformation, the A phase voltage is expressed as in (5). 
Va = ZzeroigSin(Wpt) + ZposigSin(wpt + Yj) + ZnegigSin(wyt + ¢;) (5) 


Where: 
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Op = 2nf,,p = 0,1,2 
f, =fo + foert 
f, = fo - foert 


0, =— tan? (r (Gre = 4) =1,2 


The introduction of sub-harmonic injection in the q-axis current induces oscillations in the distributed DG 
frequency during an islanding state, while such oscillations are not manifested in a grid connected distributed 
generation (GCDG) state. When the ROCOF surpasses the designated threshold, the inverter disengages from 
grid energization. 


2.2. Active islanding detection scheme II — (Complementary reactive power perturbation (CRPP)) 

In the second islanding identification approach, complementary reactive power perturbation (CRPP) 
is injected at PCC. These perturbations cause the frequency to move away from the threshold value during 
islanding scenario. However, there is no such frequency deviation during grid synchronized mode operation of 
DGs. The active and reactive load powers are given by (6) and (7). 


Proad = Pog + AP (6) 
Quoad = Qng + AQ (7) 
AP and AQ are the grid share components. The load — power components are expressed as in (8) and (9). 


v2 
Proad = = (8) 


1 
Quoad = 3V?pcc (= — 2nfC) (9) 


The resonant frequency and load quality factors are given by (10) and (11). 


f =) 1 
r 2nVLC 


Cc 
g=r/f (1) 


The frequency, f in the grid islanded mode is computed from (12). 


=+ ( usad ) Quoad (12) 
ie) ~ QPLoad 


During loss of mains, frequency fi be computed from (13)-(16). 


(10) 


Pag = 3 PEE ree (13) 


T T 1 
Qhoaa = 3V* pec (— = 2nfC) (14) 
fr Qne \? _ Qe 
- 2 (222.) + - fe ov 
owt Qnc \* QpG was.) Quoad 
= 4 Ga . — atx Ni (a) +4 - Sia] (6) 


Allowable frequency variation is between fiow and fhign, the NDZ is as presented in (17). 
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f f QLoad f f 
oO P high? 
Q (Ase ve a Pe < Q( 1g) ) (17) 
f flows —hoad f fhigh® 
DG 


Where: 


Increasing the injected reactive power reduces NDZ but causes more PQ issues in return. 


2.3. Active islanding detection scheme ITI — (Even harmonic perturbation (EHP)) 

In the proposed EHP method, even harmonics are intentionally introduced into q axis current. This 
causes ROCOF during the distributed generation's islanding state. During grid-tied operation, the grid serves 
as an effectively infinite sink, absorbing and neutralizing all injected disturbances. The generation of even 
harmonic perturbations is achieved by the application of Park and inverse Park's transformations given by (18) 
and (19). 


Cos@ ~=Cos (0 = =) Cos (0 + =) 


ud ‘ 3 ua 

| ==|-Sind -Sin(@@-=) -sin((e +=) | (18) 
ud 1 1 1 uc 

2 2 2 
2 Sin8 Cos 1 a 
F 27 27. 

| =2]sin@-) Cos(@-F) 1 : 1 (19) 
uc Sin(0 + =) Cos(@ + =) 1| luo 


3. RESPONSE OF THE PASSIVE METHOD 

After the perturbations are fed from the DG inverter, the response is to be monitored to identify the 
islanding scenario. The perturbations, used in the three active methods in this work are, injected through q axis 
current. Q axis current is a major deciding factor of reactive power from inverter. In the event of a mains failure, 
frequency change is determined by reactive power. In order to detect the islanding situation in addition to each 
of the current approaches, ROCOF is the preferred answer. | Hz/s is the threshold set for this method. Voltage 
THD is also measured and its upper limit is set at 5%. When both ROCOF and THD go above the threshold 
limits, islanding is confirmed. 


3.1. ROCOF measurement 

ROCOF based islanding detection works on the principle that the PCC frequency variation is 
manifested during the transition to an islanding state. The frequency measurement is done through a phase 
locked loop module. Rate of change of frequency calculated is passed through a low pass filter to filter out the 
noises. Then its absolute value is found out. The implementation of ROCOF is depicted in Figure 2. A trip 
signal is initiated if the measured value of ROCOF is greater than a threshold, which is | Hz/s. 


DGVoltage —> PLL > did —> UF => (|u| — ROCOF 


Figure 2. Block diagram of ROCOF measurement 


3.2. THD measurement of DG voltage 

Islanding detection methods based on harmonic distortion analysis involve the continuous monitoring 
of harmonic variations resulting from load changes in isolated DGs. Specifically, in DGs featuring inverter- 
based generation units, the operation of these inverters introduces distinct harmonic components. Detecting 
islanding scenarios can be accomplished through the measurement of harmonic distortions in either voltage or 
current signals. A viable approach for leveraging harmonic distortions in islanding detection is the computation 
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of total harmonic distortion (THD) and the monitoring of its fluctuations over time. The measurement of THD 
in DG voltage is conducted using (20), in which, V; represents the root mean square (RMS) value of the 
fundamental voltage component, while V;, signifies the RMS value of the harmonic component with index ‘h.' 
The variable 'N' denotes the maximum number of considered harmonic components. 


sh=N v2 
Ven = 1005», (20) 


Vi 


4. RESULTS AND VALIDATION 

The MATLAB simulation in Simulink uses a 50 kW DG powered by solar energy and is depicted in 
Figure 3. The output voltage of DC/DC converter is 900 V. The DG has a local parallel RLC load and is tied 
with grid at PCC. The islanding is simulated at 4.2s. The simulation details are illustrated in Table 2. The 
variation in irradiation and bus voltage are presented in Figure 4 and Figure 5. The voltages are shown in the 
Figure 6. The corresponding currents are depicted in Figure 7. 


PWM Generator 


vp) (DC-DC) ISLANDING DETECTION USING SHP, CRPP, EHP 
1 < duty IN PV INTEGRATED GRID SYSTEMS 
Pret Bigna 1 <V_PV> Goto 
Signal Builder a 
Bus Saturation 
Selector 1e-05 s. 
Powergui 


3 ae l a | From 


Temparature Py Array 


> Vbus 
IGBT/Diode 
— 7] 
Voltage Measurement Three-Phase Grid 


INVERTER Series RLC Branch 


. 


Figure 3. Simulation diagram of the IDM in grid synchronized DG 


Table 2. Test system details 


Features Ratings 

Load 50 kw 
Q factor 2.5 
f, 5,000 Hz 
Vp 245 V 
DC bus Voltage 900 V 
fo 50 

fsup 20 Hz 
cai Ta sup 1% of id-rated 
CRPP iqcrre ISA 
EHP iqeup 12A 


Resistance, R 3.4Q 
Load details Inductance, L 0.00431 H 
Capacitance,C 0.0023 F 


4.1. Sub harmonic perturbation (SHP) 

The SHP is incorporated into the q-axis current of the inverter system. Specifically, a 20 Hz SHP of 
sufficient magnitude, equal to 1% of the rated Id, is introduced into the q-axis current of the inverter as 
presented in Figure 8. The ROCOF during islanding is depicted in Figure 9. As soon as the loss of mains 
happens at 4.2s, ROCOF exceeds the threshold of 1Hz/sec, at t = 4.42s. Thereby islanding state is confirmed. 
The THD has gone above the threshold of 5%. THD of DG voltage during grid synchronized as well as 
islanding state are shown in the Figure 10(a) and Figure 10(b). The SHI induces low-frequency oscillations in 
the DG absolute frequency, when it operates in an islanded state. Notably, these oscillations are absent in grid 
connected distributed generation (GCDG) mode, as sub-harmonics are efficiently channeled directly into the 
utility grid. But in the islanded mode these oscillations are present. This observation is illustrated in 
Figure 11(a) and Figure 11(b). 
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Figure 6. DG-voltage, grid voltage, and load voltage 
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Figure 4. Variation of irradiation 
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Figure 8. SHP injected from DG inverter 
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Figure 7. DG current, grid current, and load current 
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Figure 9. ROCOF in SHP method 
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10. THD of DG voltage in SHP during (a) grid synchronized mode and (b) during islanding mode 
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Figure 11. Oscillations in the frequency waveform during islanding state (a) with SHP (b) without SHP 
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4.2. Complementary reactive power perturbation (CRPP) 
In the Figure 12, CRPP injected from the DG inverter is presented. The ROCOF is shown 
in Figure 13. The THD has beyond the 5% cutoff point as soon as the loss of mains happens at 4.2 seconds. 
ROCOF exceeds the threshold of 1Hz/s at t=4.3 seconds after the occurrence of islanding at 4.2 s. The THDs 
of DG voltage during grid synchronized as well as islanding state are shown in the Figures 14(a) and 14(b). 
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Figure 12. Complementary reactive power perturbation 
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Figure 13. The ROCOF in CRPP method 
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Figure 14. THD of DG voltage in CRPP during (a) grid synchronized mode, and (b) islanding mode 
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Even harmonics injected along with the q axis current is depicted in Figure 15. The injection of even 
harmonics causes the ROCOF to go beyond the threshold of 1Hz/s at t=4.57seconds after the occurrence of 
islanding at 4.2 s, which is shown in Figure 16. THDs of DG voltage during grid synchronized as well as 
islanding state are shown in the Figures 17(a) and 17(b). 
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Figure 15. Injected second order harmonics 
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Figure 16. The ROCOF in EHP method 
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5. COMPARISON OF SHP, CRPP AND EHP 

The comparative analysis of NDZ and FDZ of the three approaches are presented in Table 3. The 
THD values of DG voltage during grid synchronized mode for different load quality factors is presented in 
Table 4. It is also shown in Figure 18. CRPP method shows least THD during grid synchronized mode. The 
THD values of DG voltage during islanded mode for different load quality factors is presented in Table 5. It is 
also shown in Figure 19. CRPP method shows least THD during grid synchronized mode. The islanding 
detection time of three methods for various load quality factors are presented in Table 6 and is also shown in 
Figure 20. CRPP method shows least THD during both islanding mode and Grid Integrated mode. The detection 
time of the three approaches (SHP, CRPP, EHP) are compared with the existing methods and is presented in 
Figure 21. CRPP approach has the least detection time as it is closest to the center of the radar web. 


Table 3. Comparison of NDZ and FDZ 


Method NDZ FDZ 
AP=0, AQ=0 Load QF 3 phase Induction 3 phase capacitor bank 3-Phase and Non-Linear Load 
Motor Load switching switching(25k VAR) LLG faults (Z;= (Bridge Rectifier with 
(20hp) 0.01) 25kW Load) 
SHP As AP <5% NoNDZ No False tripping No False tripping No False tripping No False tripping 
Detection time goes _—_ up to 
above 2s QF=5 
CRPP No NDZ NoNDZ No False tripping No False tripping No False tripping —_No False tripping 
up to 
QF=5 
EHP No NDZ NoNDZ No False tripping No False tripping No False tripping No False tripping 
up to 
QF=5 
Table 4. THD of DG voltage in GC mode Table 5. THD % in Islanded mode 
QF=1 QF=1.5 QF=2 QF=2.5 QF=3 QF=4 QF=5 IDMs QF=0.5 QF=1 QF=1.5 QF=2 QF=2.5 QF=3 QF=4 QF=5 
0.1 0.08 0.05 0.03 0.02 0.01 0.009 EHP 12 112 109 99 9.78 9.3 9 8.5 
0.07 0.05 0.04 0.02 0.015 0.01 0.009 SHP 7.2 6.5 5.8 5.2 5.03 49 47 4,55 
0.07 0.04 0.03 —0.01 0.01 0.009 0.008 CRPP 7 6.2 555: 5.1 5.01 48 41 39 


Table 6. Islanding detection time (IDT) 


IDMs QF=0.5 QF=1 QF=1.5 QF=2 QF=2.5 QF=3 QF=4 QF=5 
EHP & 
ROCOF 0.31 0.32 0.33 0.35 0.37 0.39 04 0.42 
SHP & 
ROCOF 0.15 0.16 0.18 0.2 0.22 0.24 0.25 0.27 
CRPP & 
ROCOF 0.082 0.086 0.09 0.095 0.1 0.13 0.14 0.16 
———CRPP comm SHP & 
ROCOF 
— ——BRPP& 
QF=1 “rrr SHP ROCOF 
mmr EHP & 
\ ROCOF 
= OF=1.5 ) QF=15 
OF =2 
0F925 
Figure 18. THD% in grid Figure 19. THD% inislanded _ Figure 20. Islanding detection time 
synchronized mode mode for different load quality factors 
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Figure 21. Comparison of Islanding detection time of the three approaches with the existing methods 


CONCLUSION 
The performance assessment of three q axis perturbation based active islanding detection methods 


(SHP, CRPP, EHP) are presented in this paper. The NDZ and FDZ analysis are also presented in this paper. In 
all three approaches, ROCOF has exceeded the 1Hz/sec threshold limit when an islanding condition occurs. 
From the results it has been identified that the CRPP method is able to detect the islanding scenario with least 
detection time of 0.1seconds (QF=2.5). THD of DG voltage is least for CRPP method during grid synchronized 
mode. But this superior method is able to drift the %THD of DG voltage well beyond the threshold value of 
5%, which is a desirable feature of any active method. Therefore, the CRPP method outperforms the SHI and 
EHP methods for various load quality factors. 
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